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Abstract 

It has long been known that the multivariate integration problem for the unit ball 
in C r ([0, l] d ) is intractable for fixed finite r. H. Wozniakowski has recently conjectured 
that this is true even if r = oo. This paper establishes a partial result in this direction. 
We prove that the multivariate integration problem, for infinitely differential functions 
all of whose variables are bounded by one, is not strongly tractable. 

1 Introduction 

Multivariate integration is a classical problem of numerical analysis that has been studied for 
various normed spaces Fd of functions of d variables. In practical applications d is frequently 
very large, even in the hundreds or thousands. 

Tractability and strong tractability of multivariate integration has been recently thor- 
oughly analyzed. These concepts are defined as follows. We consider the worst case setting 
and define n(e, Bd) as the minimal number of function values that are needed to approxi- 
mate the integral of any / from the unit ball of Fd with an error threshold of e. We want to 
know how n(e, Bo) depends on d and e~ l . The problem is tractable if n(e, Bo) is bounded by 
a polynomial in e -1 and d, and strongly tractable if the bound only depends polynomially 
on e -1 , with no dependence on d. 
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The dependence of n(e, B d ) on e 1 has been studied for many years, and bounds for 
n(e, B d ) in terms of e _1 are known for many B d . For instance, let F d = C r ([0, l] d ) be the 
space of r times continuously differentiable functions defined on the ci-dimensional unit cube 
with the norm given as the supremum of the absolute values of all partial derivatives up to 
the rth order. Bakhvalov proved in 1959, see |NW01j . that there exist two positive numbers 
c r)d and C r>d such that 

c r , d e~ d/r < n(e,B d ) < C r4 e- d ' r Ve e (0,1). 

When r is fixed this implies that multivariate integration in C r ([0, l] d ) is intractable. Indeed, 
n(e, B d ) cannot possibly be bounded by a polynomial in e _1 and d since the exponent of e^ 1 
goes to infinity with d, and n(e, B d ) is exponential in d. However if r varies, that is, when we 
consider the spaces F d = C r ^ d \[0, l] d ) with sup d d/r(d) < oo, the behavior of n(e, B d ) is not 
known since we do not have sharp bounds on c T ^ d and C T ^ d . In fact, the best known bounds on 
c T)d are exponentially small in d, while the best known bounds on C r ^ d are exponentially large 
in d, see again |NW01] . Thus we can neither claim nor deny tractability or strong tractability 
in the class C r( - d \[0, l] d ) with sup d d/r(d) < oo on the basis of Bakhvalov's result. 

The conjecture formulated in [W03] states that multivariate integration in C r ^ d \[0, l] d ) 
is intractable even if r(d) = oo. That is, even when we consider infinitely differentiable 
functions with all partial derivatives bounded by one, n(e, B d ) cannot be bounded by a 
polynomial in e" 1 and d. Although we are not able to establish this conjecture in full 
generality, we shall prove that multivariate integration in C°°([0, l] d ) is not strongly tractable. 
This is achieved by showing that for a fixed n, the nth minimal error goes to one as d 
approaches infinity. More precisely, we show that for any n and 77 there exists d = d(n, 77) 
such that for any linear algorithm there exists a polynomial that is a sum of univariate 
polynomials, which belongs to the unit ball of C°°([0, l] d ), whose integral is at least 1 — T], 
and the algorithm outputs zero. This proof technique allows us to prove the lack of strong 
tractability, but seems too weak to establish the lack of tractability. 

2 Multivariate Integration in C°°([0, l] d ) 

We precisely define the problem of multivariate integration studied in this paper. Let F d = 
C°°([0, l] d ) be the space of real functions defined on the unit cube [0, l] d that are infinitely 
differentiable with the norm 

\\f\\ d = sup{ \D a f(x)\ : x G [0,l] d , a - any multiindex }. 
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Here a = [«i, a 2 , ■ ■ ■ , a d ] with non-negative integers atj, \a\ = a>i + a 2 + ■ ■ ■ + a d , and 

£> a f = f 



dx^ 1 2 . . . <9x 



ay • 

(2 



stands for the differentiation operator. 

Let B d denote the unit ball of this space. The multivariate integration problem is defined 
as an approximation of integrals 



W)= [ f(t)dt WfeB d , 

J\0,l] d 



'[0,1] 

where the integral is taken with respect to the Lebesgue measure. 

It is well known that adaption and the use of non-linear algorithms do not help for 
the multivariate integration problem, as proven in |B71j . see also [NWOlj . That is why we 
consider only linear algorithms, 

n 

An,d(f) = ^2 a jf( X j) 

i=i 

for some real coefficients aj and some sample points Xj from [0, l] d . Here n denotes the 
number of function values used by the algorithm. Of course the ctj and Xi may depend on n. 
The (worst case) error of the algorithm A n ^ is defined as 

err(A n4 ) = sup{|/ d (/) - A n4 (f)\ : / G B d }. 

Let LIN n <d denote the class of all linear algorithms that use n function values. The nth 
minimal error is defined as 

e(n,B d ) = inf{err(v4 nid ) : A n4 e LIN njd }. 

We shall prove the following theorem: 

Theorem 2.1 For any positive integer n we have 

Jim e(n,C°°([0,l] d )) = l. 

0— >oo 

This theorem easily implies that multivariate integration in C°°([0, l] d ) is not strongly 
tractable. Indeed, were it strongly tractable, we would have a polynomial bound on n(e, B d ) 
independent of d, thus having a linear algorithm of error at most e and using at most 
n = n(e, B d ) function values. Taking, say, £ < \ and n > Bd) we would get e(n, B d ) < | 
independently of d, and this would contradict the theorem that e(n, B d ) goes to one when d 
approaches infinity. 
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2.1 Proof of the Theorem 



We take an arbitrary positive integer n and rj G (0, 1). The idea of the proof is to separate 
variables and, for sufficiently large d and any A n ^ G LIN ni( j, to find a polynomial / G Bd 
that is a sum of univariate polynomials such that \Id{f) ~ Ai,d(/)| > 1 — It will suffice to 
find such a polynomial for which f(xj) = at all points Xj used by A nd . Then A Ujd (f) = 0, 
and thus err(A n ^) > > 1 — V- Since ^ is an arbitrary linear algorithm this implies 

that e(n, Bd) > 1 — r\ for sufficiently large d. For the zero algorithm A U)d = we have we 
have err(y4 nj d) = 1 and hence we have e(n,Bd) < 1. Since 77 can be arbitrarily small, this 
completes the proof of Theorem 12.11 

Suppose for the moment that we have the following lemma, which will be proven in the 
next subsection. 

Lemma 2.2 For any positive integer n, and any rj G (0, 1) there exists a constant K v n such 
that for any choice ofyi,yz, . . . ,y n G [0, 1] there exists a polynomial f : [0, 1]— >R satisfying 
the following conditions: 

1. max^o,!] |/(ar)| < 1, 

2. maxjfe =0) i... maxa-gici] \ f {k) {x)\ < K VtU , 

3. j'l f(x) dx>\-r], 

4. f( Vj ) = for j = 1,2,..., n. 

Having Lemma [2.21 we take any d > K v ^ n and any A G LIN nj d that uses sample points 

" dl G [0,l] d 



[rf / >"» *^ ■ • ' ' 

3 ' 3 ' ' ' ' ' i J 



for j = 1,2,..., n. 

For i = 1,2, ... ,d, let f\ be the polynomial given by Lemma 12.21 for yj = x*, with 
j = 1, 2, . . . , n. Consider the multivariate polynomial 



1 d 

f{t u t 2 ,...,t d ) = -j^lfi(U) U G [0,1]. 

1=1 



The values of / are bounded by 1 since they are arithmetic means of the values of f\ from 
[—1,1]. Any mixed derivative of such a function / is 0, while 



Qk 

dx) 



d 



f 



< 



K 



< 1. 
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Thus / belongs to B d . Additionally, 



I f(t)dt = -J2[ fi(x) dx > 1 - r). 
J[o,i] d a tX J l A 

Furthermore f(xj) = since fi(x^) = for all j = 1,2, ...,n. Thus / is a function we 
needed to prove Theorem 1 2. 1[ 

2.2 Proof of the Lemma 

We will use the Stone- Weierstrass theorem to find a function satisfying Lemma 12.21 

For 5 = r//(7n), let the function g : [0, 2+5] ->R be defined as 1-5 on [0, l-S\U[l+S, 2+5], 
—25 at 1 and linear on [1 — 5, 1] and [1, 1 + 5]. It is obviously a continuous function, so by the 
Stone- Weierstrass theorem we can approximate it by a polynomial P of degree N = N(rj, n) 
such that 

max \g(x) — P(x)\ < 5. 

The polynomial P is negative at 1 and positive at 1 + 5, so it has a root at some yo G 
(1, 1 + 5). Let Pi(x) = P(x + yo — Hi). As y — Ui G (0, 1 + 5) the polynomial Pj satisfies 
max xg [ 0il ] \Pi(x) - g(x + y - y-i)\ < 5. Now take 

n 

f(x) = H P ^ Va;G[0,l]. 
i=i 

Note that f(yj) = Yi7=i ^ 2 (yj + 2/o — = since the jth factor is P 2 (yo) = 0. 
The polynomial Pi satisfies 

1 - 25 < P&) < 1 Va; G [0, 1] \ [ Vi - 25, Vi + 6). 

Thus 

n 

(1 - 25) 2n < f(x) < 1 Va; G [0, 1] \ \J[ yi - 25, y t + 5] 

8=1 

and, of course, f(x) > on the whole interval [0,1]. This allows us to approximate the 
integral of /. Indeed, 

/ f(x)dx> [ f{x)dx > (l-35n)(l-25) 2n . 
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Using the Bernoulli inequality we conclude that the last expression is at least 

(1 - 3rwJ)(l - 4n<5) > 1 - 7n<J = 1 - 77. 

The function / is a polynomial on [0,1]. Its coefficients are continuous functions of 
(yi, y2, ... , y n ) £ [0, l] n since the coefficients of each are continuous functions of y^ and / 
is the product of Pj's. The upper bound of the jth derivative of / is a continuous function of 
/'s coefficients, and thus a continuous function of (yi,yz, ■ ■ ■ ,y n )- As a continuous function 
on a compact set it is bounded for each j, and so all derivarives up to the 2niVth order have 
a common bound, say, K v>n , independent of (y^ y 2 , ■ ■ ■ , y n ). This means that the second 
condition of Lemma 12.21 is satisfied and the proof of Lemma 12.21 is completed. 
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